Spin-glass behavior in the random-anisotropy Heisenberg model 
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We perform Monte Carlo simulations in a random anisotropy magnet at a intermediate exchange 
to anisotropy ratio. We focus on the out of equilibrium relaxation after a sudden quenching in the 
low temperature phase, well below the freezing one. By analyzing both the aging dynamics and 
the violation of the Fluctuation Dissipation relation we found strong evidence of a spin-glass like 
behavior. In fact, our results are qualitatively similar to those experimentally obtained recently in 
a Heisenberg-like real spin glass. 
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INTRODUCTION 

Random magnetic anisotropy seems to be a fundamen- 
tal ingredient for any realistic description of amorphous 
materials, which are systems of both practical and the- 
orical relevancy. For the particular case of amorphous 
alloys 0,0 of non-S-state rare earths and transition met- 
als (RE-TM), such as Tb^Fei-^, the three dimensional 
classical Heisenberg model with random uniaxial single- 
site anisotropy (RAM) is considered to be the proper 
model for studying both their thermodynamic and dy- 
namic properties. This model was introduced by Harris, 
Plischke, and Zuckermann, @] who performed a mean- 
field calculation and found a ferromagnetic (FM) phase 
at low temperature. Later on, Pelcovits, Pytte and Rud- 
nick Q claimed, using an argument similar to that used 
by Imry and Ma 0] for the random-field case, that such a 
FM phase is not stable in three dimensional RAM model, 
for any finite value of the anisotropy. Since then, the na- 
ture of the ordered phase at low temperatures and its 
dependence on the degree of anisotropy is a subject of 
controversy. Recently, Itakura || proposed, by using 
Monte Carlo simulations and referring to former works in 
the literature, a schematic phase diagram for the RAM 
where different kinds of order can be found, depending 
both on the temperature and the anisotropy strength of 
the system. 

When an amorphous material is cooled, it can even- 
tually get blocked at certain temperature T/, at which 
the magnetic moments freeze pointing in random direc- 
tions. The value of Tf strongly depends on the degree of 
anisotropy, the strength of the interactions between do- 
mains and the cooling rate. It is worth here to stress 
that this freezing process is a dynamical phenomenon 
which can not be associated to any true thermodynam- 
ical phase transition. In particular, this phenomenology 
has been reported, during the last years, in the study of 
hard magnetic amorphous alloys 0,0; which has been 
also successfully simulated using a slightly modified ver- 
sion of the RAM Q. On the other hand, in experimental 
spin glasses 0], zero field cooling (ZFC) and field cooling 



(FC) curves of magnetization versus temperature are use- 
ful to estimate the characteristic freezing temperatures of 
the systems. 



Most of the numerical effort in the study of the RAM 
model concerns the case of strong or infinite anisotropy. 
In this case the model seems to present a low temperature 
spin glass like phase (usually called speromagnetic ^H). 
On the other hand, in the weak anisotropy limit, the sys- 
tem tends to order locally in a ferromagnetic state |l2( 
(also called asperomagnetic) . When increased to inter- 
mediate values, the anisotropy seems to destroy the asper- 
omagnetic ferromagnetic state, and the systems orders in 
the so called correlated speromagnetic or correlated spin- 
glass. It has been recently verified in Ref.@, by means of 
a very extensive numerical simulation, that the compe- 
tition between exchange and anisotropy gives place to a 
quasi-long-range order (QLRO) low temperature phase 
characterized by frozen power-law spin-spin correlations. 
It is important here to remark that all the theoretical re- 
sults ultimately confirm the observed lack of long range 
ferromagnetic order observed experimentally in magnetic 
materials with isotropic easy axis distribution [lij . 



Although the RAM model has deserve much attention 
during the last years, most of the works were concerned 
on its equilibrium properties as well as on its magnetic 
behavior, paying little attention to its relaxation dynam- 
ics. The main question we want to address in this work 
concerns the possible existence of a spin-glass like dy- 
namical behavior associated with QLRO low temperature 
phase in the intermediate anisotropy regime. We analyze, 
through Monte Carlo simulations, the out of equilibrium 
relaxation of the three dimensional RAM model defined 
on a cubic lattice. In particular, we report results ob- 
tained in the low temperature phase (well below Tf) and 
for intermediate values of the anisotropy to exchange ra- 
tio. 
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THE MODEL 

The system is ruled by the following classical Heisen- 
berg Hamiltonian: 

H=-J ^ Si- S, -Dj2(fii- - ^ ■ Si (1) 

<i,j> i 

where D and J > are the anisotropy and the exchange 
strength, respectively, and Hi is an external field act- 
ing on the site i. The spin variable Si is a three com- 
ponents unit vector associated to the i-th node of the 
lattice and the first sum runs over all nearest-neighbor 
pairs of spins. Hi is unit random vector that defines the 
local easy axis direction of the anisotropy at site i. These 
easy axis are quenched variables chosen from a isotropic 
distribution on the unit sphere. The simulations were 
performed in a system of N = L 3 spins (L = 15), using 
a Monte Carlo Metropolis algorithm |14j with periodic 
boundary conditions (in fact, we have compare different 
system sizes up to L = 20, confirming that for L = 15 
finite size effects become very small). The ratio between 
the anisotropy and the exchange strength was fixed at 
the value D / J = 3.5, which is comparable to those values 
observed in real amorphous material and clusters ferro- 
magnetic alloys Following the ideas used in Ref. 0, 
at each spin actualization the direction of the spin is ad- 
justed in such a way to maintain an acceptation rate close 
to 0.46. 



RESULTS 

The first step was to locate the freezing temperature Tf 
for these particular values of the parameters by looking 
for the temperature at which zero field cooling (ZFC) and 
field cooling (FC) curves split each other. FigQ] shows 
the magnetization divided by the field M/H (dc suscep- 
tibility) as a function of T j J, both for the ZFC and FC 
processes. The simulation protocol used is the following: 
we performed 1000 Monte Carlo steps (MCS) at a given 
temperature and at a constant acceptance rate in order 
to equilibrate the system, after which, we used other 1000 
MCS to get a time average of the magnetization before 
decreasing the temperature. In all the simulations pre- 
sented in this work we have used between 20 and 40 dif- 
ferent realizations of the disorder to average the results. 
We find that the freezing temperature Tf/J is close to 
0.5. At very low temperature the dc susceptibility in the 
ZFC curves is constant and the magnetization almost 
zero, indicating a speromagnetic spin-glass like order. In 
other words, the spins are frozen into random orienta- 
tions with average correlation over at most one lattice 
parameter (we have verified this behavior by calculating 
the dependence of the spin-spin correlation function as 
function of the distance). As temperature increases the 
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FIG. 1: Magnetization divided by the field M/H (dc sus- 
ceptibility) as a function of the rescaled temperature T / J for 
different applied fields, both for the ZFC and FC processes. 

dc susceptibility also increases, indicating that the sys- 
tem goes to an anisotropic asperomagnetic phase. Here 
again, the same conclusion has been verified by analyz- 
ing the correlation function which stabilizes in a non zero 
value for distances larger than approximately four lattice 
units. 

On the other hand, it can be seen that as T — > 0, 
Mfc — * 0.5, as expected in an asperomagnetic phase 
(and in agreement with other theoretical |l2| and nu- 
merical predictions ^3|). Due to local character of this 
asperomagnetic phase, a sufficiently large system could 
eventually develop a local order without any global mag- 
netization (QLRO), as in fact is observed in amorphous 
materials. Notice that both curves present, above Tf, an 
inflection point which is very close to the Heisenberg crit- 
ical temperature, indicating the loss of magnetic order 
and the entrance in the paramagnetic phase (or super- 
paramagnetic phase, in particulate systems). It is worth 
here to remark that the FC curves display a clear cusp, 
that closely resembles the experimental results obtained 
in spin glasses, as for intance in AgMn confirming 
the random freezing of the spin orientations. 

Disordered systems, when suddenly quenched down 
into the low temperature phase, suffer a drastic slowing 
down of their relaxation dynamics. At the same time, 
a very strong dependence on the history of the sample 
emerges, a phenomenon usually called aging. In a real 
experiment, the simplest way for measuring aging is by 
suddenly quenching the system without field into the or- 
dered phase. The system ages in this phase during certain 
waiting time t w , at which the field is switched off. The 
measurement of the relaxation of the magnetization then 
strongly depends on both, the age or waiting time t w and 
the time t elapsed since the field was turned off, indicat- 
ing the loss of time translation invariance (TTI) proper 
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FIG. 2: Correlation curves C(t w + t,t w ) as function of time 
t (MCS) obtained at T/J = 0.2, for L = 15 and for different 
values of t w . 



of any equilibrated state. In a computer simulation the 
same effect can be visualized by measuring the two time 
auto correlation function after a sudden quench from in- 
finite temperature into the low temperature phase: 



N ^ 



Si(t w ) ■ Si(t + t w ) 



(2) 



where t w is the time elapsed after the quenching and [• • •] 
means an average over the disorder. 

In Fig. we present the curves of C(t w + t,t w ) ob- 
tained at T/J = 0.2, L = 15 and for different values of 
t w . The plot clearly confirms the appearance of aging, 
characterized both by the loss of TTI and the fact that 
the system decays slower as its age t w increases. 

Actually, aging is so ubiquitous in nature, that one can 
wander whether it is useful or not to look for this phe- 
nomenology. But fortunately, the peculiar dependence of 
C(t w +t,t w ) on t and t w on a great variety of systems 
(both theoretical and real systems) suggests the existence 
of only a few universality classes associated to the out of 
equilibrium relaxation of the model, as occurs, for in- 
stance, in critical phenomena and coarsening dynamics. 
All this indicates that, despite any microscopic difference 
between different systems, the relaxation must be domi- 
nated by a few relevant ingredients. 

In order to classify the universality class, we looked 
for the best data collapse of the curves obtained for 
different waiting times, checking different scaling forms. 
The best results we obtained are presented in Figure 
where we followed the standard scaling procedure used 
in spin glasses materials (like for instance CdIno.3Cr1.7S4 
[lflj)- The two time autocorrelation function is assumed 
to have two different components C(t w +t, t w ) = C s t{t) + 
C a g{t w +t,t w ), the first one corresponding to an station- 
ary part (independent of t w ) and the second part (the 



FIG. 3: Rescaled correlations curves for three different values 
of T/J (0.1, 0.15 and 0.2). Each curve we shown in the graph 
correspond to four correlation curves C(t w +t,t w ) with t w = 
5000, 10000, 20000 and 40000 (measured in MCS). 



TABLE I: Fitting parameters of the scaling. 
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aging one) depending both on t and t w . The stationary 
part is a power law function A(r /t) a , which is predomi- 
nant at very small times. The aging part is a function of 

A(i, t w )/t w , where 



\ = t w [(l+t/t w ) 1 -> i -l}/[l- tJ ] 



(3) 



The curves presented in FigO correspond to three dif- 
ferent values of T/J (0.1, 0.15 and 0.2), all of them well 
below the freezing temperature of the system. In each 
curve we have collapsed the four different curves obtained 
for four different waiting times t w = 5000, 10000, 20000 
and 40000 (measured in MCS). Note the excellent super- 
positions obtained, which extends over the complete time 
span simulated. It is worth to stress that this scaling law 
has been vastly used in the study of aging dynamics in 
real spin glass materials, yielding excellent data collapses 
of the experimental results 

The set of fitting parameters obtained (shown in table 
HJ| are in a good agreement with those obtained experi- 
mentally in real spin glasses ^^|. This behavior indicates 
that C(t w + 1, t w ) scales as t/t^ at short t (since fi < 1), 
while in the large t limit and for the values of \i obtained 
(which are all close to 1) its behaviors is almost logarith- 
mic on t/t w (as expected in a activated scenario 



igaritJ 

0). 
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Finally, we analyzed the Fluctuation Dissipation Rela- 
tion (FDR), which can be expressed as |2lj |: 



X(t w + t, t w ) dC(t w + t, t w ) 



3T 



8ty 



(4) 



where R{t w + 1, t w ) = l/N £. d (Sf(t w + t)) /dhi(t w ) 
is the response to an external magnetic field hi(t) in the 
z direction and X(t w + t,t w ) is the fluctuation dissi- 
pation factor. In equilibrium the Fluctuation Dissipa- 
tion Theorem (FDT) holds and T X(t w + t, t w ) = 1, 
while out of equilibrium X depends on t and t w in 
an non trivial way. It has been conjectured HH that 
X(t w + t,t w ) = X[C(t w + t,t w )]. This conjecture has 
proved valid in all systems studied to date. 

Instead of considering the response function it is easier 
to analyze the integrated response function 



t m +t 



R{t w + t, s) ds. 



(5) 



Assuming X(t w + t, t w ) = X[C(t w + t, t w )\ one obtains 
3T X {t w +t,t w )= f X{C)dC (6) 

JC{t w +t,t w ) 

and by plotting 3Tx vs. C one can extract X from the 
slope of the curve. In particular, if the FDT holds X = 1 
and 3 Tx(t) = (1— C(i)); any departure from this straight 
line brings information about the non-equilibrium pro- 
cess. In numerical simulations of spin glass an d 
structural glass models [13, S| it has been found that, 
in the non-equilibrium regime, this curve follows another 
straight line with smaller (in absolute value) slope when 
t/t w ^> 1. In this case the FD factor X can be interpreted 
as an effective inverse temperature (2^| T e // —T/X. At 
time t w we took a copy of the system spin configuration, 
to which a random magnetic field hi(t) = hu was ap- 
plied, in order to avoid favoring the QLRO [2g, |2j|; a 
was taken from a bimodal distribution (e; = ±1). Us- 
ing the results from the FC and ZFC calculations, the 
strength h of the field was taken small enough to guaran- 
tee linear response; the integrated response then equals 
x(t w +t,t w ) = m(t w +t,t w )/h, where m(t w +t,t w ) is 
the staggered magnetization conjugated to the field hi(t), 
averaged over the random field variables. 

In Fig. Elwe display Tm(t w +t, t w )/h vs. C(t w +t, t w ) 
in a parametric plot. The curves correspond to the three 
different values of T/J plotted in Fig. in each case 
we present the results obtained for two different waiting 
times t w = 5000 and 40000 MCS. In all the cases studied 
we observe a typical two time scale separation behav- 
ior, proper of real spin glasses. At t = the system 
starts in the right bottom corner (fully correlated and 
demagnetized) and during certain time (that depends in 
this case both on t w and /i) it runs over the equilibrium 



0.08 



0.06 



0.04 



0.02 



0.00 



T/J=0.2 \ 








" - T/J=0. 15' - 






~ ■-'**■« \ 






t w = 5000 MCS 




t™ = 40000 MCS 


T/j=o*r^\ 



0.5 



0.6 



0.7 0.8 

C 



0.9 



1.0 



FIG. 4: Parametric plot of Tm(t + t w ,t w )/h versus C(t + 
t w ,t w ) for three reduced temperatures T/J (0.1, 0.15 and 
0.2 from bottom to top) and two different waiting times t w 
(5000 and 40000 MCS). The solid line indicates the validity 
of Fluctuation Dissipation Theorem (slope 1/3). 



straight line, indicating the existence of a quasi equilib- 
rium regime. Nevertheless, at certain time the system 
clearly departures from this quasi-equilibrium curve and 
moves along a different straight line, but with a different 
(smaller) slope, indicating an effective temperature that 
is larger than the temperature of the thermal bath. This 
one-step temporal regime observed in this model is very 
common in structural glasses but also in spin glass ma- 
terials. Both numerical results on the Heisenberg spin 
glass with weak anisotropy [28] and experimental mea- 
surements in CdIno.3Cr1.7S4 spin glasses, present this 
kind of dynamical behavior. 



CONCLUSIONS 

In this paper we have studied the out of equilibrium 
dynamics of the RAM model. The parameters have been 
chosen in such a way to model systems with a freezing 
temperature well below the ordering temperature. We 
have restricted ourselves to consider the case of interme- 
diate values oiD/J, where both effects exchange and ran- 
domness, compete with each other. This peculiar regime 
is specially interesting, since on one hand there exists 
certain degree of controversy about the expected order- 
ing of the system, and on the other hand it allows to 
model different interesting real magnetic materials. Our 
analysis, based the study of ZFC-FC curves, aging and 
on the FDR in the low temperature phase, are consistent 
with the existence of a spin-glass ground state of the 
model, where the slowing mechanisms are then related 
to the topology of the energy landscape of the model. 
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